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Abstract 


This  report  describes  the  test  analysis  software  developed  at  Carleton  University  (CEW)  that  was 
written  to  support  the  work  on  the  Conducted  Energy  Weapons  (CEW)  Test  Procedure,  version 
1.1  (published  31  July  2010).  This  report  documents  the  software,  its  structure,  its  business  rules, 
and  it  identifies  any  data  elements  required  to  fully  support  the  published  test  procedure.  This 
report  is  designed  to  support  the  Conducted  Energy  Weapons  Strategic  Initiative  (CEWSI),  which 
plans  to  test  a  number  of  CEWs  and  do  analysis  on  legacy  test  data  with  the  objective  of  better 
understanding  and  interpreting  the  performance  of  the  devices,  with  a  goal  of  developing  an 
updated  CEW  Test  Procedure  and  providing  advice  to  Canadian  Eaw  Enforcement  on  these 
devices. 


Resume 


Ee  present  rapport  decrit  le  logiciel  d’analyse  d’essai  mis  au  point  par  I’universite  Carleton  (AI) 
qui  a  etc  elabore  pour  appuyer  les  travaux  menes  sur  les  procedures  de  deroulement  des  essais  sur 
les  armes  a  impulsions  (AI),  version  1.1  (publiees  le  31  juillet  2010).  Ee  present  rapport  decrit  le 
logiciel,  sa  structure  et  ses  regies  administratives  et  il  repertorie  egalement  tous  les  elements  de 
donnees  necessaires  pour  appuyer  entierement  les  procedures  d’essai  publiees.  Ee  present  rapport 
est  confu  pour  appuyer  TInitiative  strategique  sur  les  armes  a  impulsions  (ISAI),  qui  prevoit  la 
mise  a  Tessai  d’un  certain  nombre  d’AI  et  de  faire  Tanalyse  des  donnees  d’essai  des  anciens 
systemes  dans  le  but  de  mieux  comprendre  et  de  mieux  interpreter  le  comportement  de  ces 
dispositifs,  avec  comme  objectif  de  mettre  au  point  de  nouvelles  procedures  d’essai  mises  a  jour 
des  AI  et  de  pouvoir  conseiller  les  organismes  canadiens  charges  de  Tapplication  de  la  loi  sur  ces 
dispositifs. 
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Executive  summary 


Documentation  of  the  Carleton  University  Conducted  Energy 
Weapons  (CEW)  Test  Analysis  Software 

Andy  Adler;  Owen  Marsh;  David  P.  Dawson;  DRDC  CSS  CR  2011-18;  Defence 
R&D  Canada  -  CSS;  October  2011. 

Introduction  or  background:  There  is  a  need  for  reliable  uniform  testing  of  Conducted  Energy 
Weapons  (CEWs)  independent  of  the  manufacturer,  as  identified  by  several  studies  including  the 
Braidwood  Commission  report,  the  Report  of  the  Standing  Committee  on  Public  Safety  and 
National  Security  of  the  Conducted  Energy  Weapon.  Recognizing  this  need,  a  team  at  Carleton 
University  (lead  by  Professor  A.  Adler  and  Mr.  D.P.  Dawson)  began  a  research  program  to 
develop  techniques  and  data  analysis  methods  for  data  from  CEWs.  This  research  has  led  to: 

1)  development  of  CEW  data  acquisition  and  test  systems,  2)  Development  of  software  to 
analyse  CEW,  3)  organizing  two  Canadian  CEW  workshops,  and,  based  on  consensus  at  the 
workshop,  4)  leading  of  a  team  to  write  a  CEW  Test  Procedure  (version  1.1,  July  31,  2010). 

Results:  This  report  describes  the  test  analysis  software  developed  at  Carleton  University  (CEW) 
that  was  written  to  support  the  work  on  the  Conducted  Energy  Weapons  (CEW)  Test  Procedure. 
This  report  documents  the  software,  its  structure,  its  business  rules,  and  it  identifies  any  data 
elements  required  to  fully  support  the  published  test  procedure 

Significance:  This  report  is  designed  to  support  the  Conducted  Energy  Weapons  Strategic 
Initiative  (CEWSI),  which  plans  to  test  a  number  of  CEWs  and  do  analysis  on  legacy  test  data  with 
the  objective  of  better  understanding  and  interpreting  the  performance  of  the  devices. 

Future  plans:  CEWSI  plans  to  use  this  report  to  further  the  goal  of  developing  an  updated  CEW 
Test  Procedure  and  providing  advice  to  Canadian  Law  Enforcement  on  these  devices. 
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Sommaire 


Documentation  of  the  Carleton  University  Conducted  Energy 
Weapons  (CEW)  Test  Anaiysis  Software:  Documentation 
reiative  au  iogiciei  d’anaiyse  d’essai  eiabore  par  i’universite 
Carieton  pour  ies  armes  a  impuisions  (Ai) 

Andy  Adler;  Owen  Marsh;  David  P.  Dawson;  DRDC  CSS  CR  2011-18;  R&D 
pour  la  defense  Canada  -  CSS;  octobre  2011. 

Introduction  on  contexte:  Tel  que  cela  a  ete  revele  par  plusieurs  etudes  sur  les  armes  a 
impulsions  (AI),  notamment  le  rapport  de  la  Commission  Braidwood  et  le  rapport  du  Comite 
permanent  de  la  securite  publique  et  nationale,  11  est  necessaire  de  disposer  de  moyens 
d'evaluation  des  armes  a  impulsions  (AI)  qui  soient  uniformes  et  fiables;  ces  moyens 
d'evaluation  doivent  aussi  etre  dissocies  du  fabricant  de  ces  armes.  Conscient  de  ce  besoin,  une 
equipe  de  I'universite  Carleton  (dirigee  par  le  Professeur  A.  Adler  et  par  M.  D.P.  Dawson)  a 
entrepris  un  programme  de  recherche  visant  a  mettre  au  point  des  techniques  et  des  methodes 
d'analyse  des  donnees  pour  les  AI.  Cette  recherche  a  mene  a  :  1)  la  mise  au  point  de  systemes 
d'acquisition  de  donnees  et  d'essai  pour  les  AI,  2)  I'elaboration  d'un  Iogiciei  permettant 
d'analyser  les  AI,  3)  I'organisation  de  deux  ateliers  canadiens  sur  les  AI  et,  en  function  du 
consensus  etabli  lors  de  I'un  de  ces  ateliers,  4)  la  creation  d'une  equipe  ayant  pour  but  de 
rediger  des  procedures  de  mise  a  I'essai  des  AI  (version  1.1,  31  juillet  2010). 

Resultats:  Le  present  rapport  decrit  le  Iogiciei  de  mise  a  I'essai  congu  a  I'universite  Carleton  (AI) 
afin  d'appuyer  les  travaux  sur  les  procedures  de  mis  a  I'essai  des  armes  a  impulsions  (AI).  Le 
present  rapport  decrit  le  Iogiciei,  sa  structure  et  ses  regies  administratives;  il  repertorie 
egalement  tous  les  elements  de  donnees  necessaires  pour  appuyer  entierement  les  procedures 
de  mise  a  I'essai  publiees. 

Importance:  Le  present  rapport  est  congu  pour  appuyer  I'Initiative  strategique  sur  les  armes  a 
impulsions  (ISAI),  qui  prevoit  la  mise  a  I'essai  d'un  certain  nombre  d'AI  et  de  faire  I'analyse  des 
donnees  d'essai  des  anciens  systemes  dans  le  but  de  mieux  comprendre  et  de  mieux  interpreter 
le  comportement  de  ces  dispositifs. 

Perspectives:  'ISAI  prevoit  utiliser  le  present  rapport  pour  appuyer  I'objectif  visant  a  elaborer  et 
a  mettre  a  jour  de  nouvelles  procedures  de  mise  a  I'essai  des  AI  et  a  conseiller  les  organismes 
canadiens  charges  de  I'application  de  la  loi  relativement  a  ces  dispositifs. 
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1  Introduction 


1.1  Background 

Several  studies  including  the  Braidwood  Commission  report[l],  the  Report  of  the  Standing 
Committee  on  Public  Safety  and  National  Security  of  the  Conducted  Energy  Weapon,  the  report 
of  the  Commission  for  Public  Complaints  against  the  RCMP  and  other  provincial  reports  and 
coroners'  recommendations  have  discussed  the  need  for  reliable  uniform  testing  of  Conducted 
Energy  Weapons  (CEWs)  independent  of  the  manufacturer.  Recognizing  this  need,  a  team  at 
Carleton  University  (lead  by  Professor  A.  Adler  and  Mr.  D.P.  Dawson)  begun  a  research  program 
to  develop  techniques  and  data  analysis  methods  for  data  from  CEWs.  This  research  has  led  to 
several  outputs:  1)  development  of  CEW  data  acquisition  and  test  systems  [2],  2)  Development 
of  software  to  analyse  CEW  3)  organizing  two  Canadian  CEW  workshops,  and,  based  on 
consensus  at  the  workshop,  4)  leading  of  a  team  to  write  a  CEW  Test  Procedure  [3]. 

This  Test  Procedure  is  designed  to  enable  organizations  across  Canada  to  test  CEWs  in  a  reliable 
repeatable  manner  to  determine  whether  they  are  operating  within  manufacturer's 
specifications.  Test  results  so  obtained  will  be  usable  in  various  ways.  The  test  procedure 
contains  a  set  of  recommendations  for  measurement  of  the  performance  characteristics  of 
conducted  energy  weapons.  It  represents  the  opinions  of  its  authors,  a  group  of  subject  matter 
experts  who  have  been  involved  in  research  on  or  testing  of  CEWs,  and  is  subject  to  the 
disclaimer  (presented  in  section  0.0  of  [3]). 

To  test  CEW  weapons  data,  an  implementation  of  the  analysis  software  has  been  created  by 
Carleton  University.  This  software  may  be  used  in  the  analysis  of  the  stored  data.  It  is  available 
under  an  open  source  license  from  Dr.  Andy  Adler,  Systems  and  Computer  Engineering,  Carleton 
University. 

This  report  documents  the  Carleton  University  Conducted  Energy  Weapons  (CEW)  Test  analysis 
software  that  was  written  to  support  the  work  on  the  CEW  Test  Procedure  [3].  This  report  is 
designed  to  support  the  Conducted  Energy  Weapons  Strategic  Initiative  (CEWSI),  a  project 
implemented  within  the  Canadian  Police  Research  Centre  and  managed  by  Defence  Research 
Development  Canada  (DRDC)  through  the  Centre  for  Security  Science  (CSS).  The  CEWSI  project 
will  be  testing  a  number  of  CEWs  and  doing  analysis  on  legacy  test  data  with  the  objective  of 
better  understanding  and  interpreting  the  performance  of  the  devices,  developing  an  updated 
CEW  Test  Procedure  and  providing  advice  to  Canadian  Law  Enforcement.  A  team  of  Carleton 
University  researchers,  led  by  Professor  A.  Adler,  have  conducted  testing  on  CEWs,  and 
developed  a  software  package  to  capture,  store  and  analyze  the  data.  Additionally,  this  team 
were  lead  authors  of  the  current  CEW  Test  Procedure  [3]  as  well  as  other  papers.  In  order  to 
help  prepare  CEWSI  for  future  steps,  a  contract  for  this  work  was  left  to  Carleton  University  to 
document  the  software,  its  business  rules  and  the  data  elements  required  to  support  the  Test 
Procedure  [3]. 
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1.2  Scope 


This  document  aims  to  document  the  Carleton  University  CEW  analysis  software,  with  a  focus  on: 

1.  Development  of  a  logical  data  model  and  physical  schema  of  the  software  data  elements 

2.  Document  algorithms  and  business  rules  incorporated  into  the  software 

3.  Identify  the  additional  data  elements  and  business  rules  that  will  be  required  to  fully 
implement  the  data  requirements  of  CEW  Test  Procedure  [3], 

1.3  Conducted  Energy  Weapons  and  Electrical  Testing 

Conducted  Energy  Weapons  (CEWs)  are  increasingly  purchased  by  police  in  many  countries  as  a 
less-lethal  force  option.  The  most  widely  used  CEWs  are  the  M26  and  X26  models  from  TASER 
International  (Tl),  which  work  by  firing  two  small  darts  (electrodes)  attached  to  wires  into  the 
subject.  The  weapon  then  sends  a  pulsatile  electrical  current  into  the  subject,  designed  to  be 
sufficiently  rapid  and  with  enough  energy  to  cause  muscular  incapacitation  [4].  This  incapacitation 
is  intended  to  give  an  arresting  officer  time  to  gain  control  over  the  subject. 

Electrical  Testing  of  a  CEW  is  designed  to  measure  the  electrical  output  characteristics  of  the 
weapon  into  a  standardized  load,  and  to  calculate  parameters  to  describe  the  electrical  output 
which  are  relevant  to  the  safety  and  physiological  effectiveness  of  the  weapon.  A  schematic 
diagram  of  a  CEW  attached  to  a  test  configuration  is  shown  in  figure  1: 


Figure  1  Schematic  Diagram  of  a  CEW  attached  to  a  test  System  (reproduced  from  [3]) 
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The  design  specifications  for  the  electronics  hardware  and  analysis  software  in  [3]  are  driven  by 
a  requirement  to  limit  errors  due  to  the  data  acquisition  and  processing  to  <  1%.  However,  it  is 
assumed  that  current  high  voltage  measurement  probes  (including  voltage  and  current  probes) 
can  introduce  errors  up  to  3%.  Thus,  the  specifications  of  [3]  were  will  result  in  only  a  small 
increase  in  errors  due  to  the  front  end  measurement  probes. 

The  Carleton  University  team  has  been  using  a  CEW  measurement  system  based  on  a  voltage 
probe  (Tektronix  lOOOx  probe  6015A).  Since  December  2010  Carleton  University  has  been  using 
a  PicoScope  dual  channel  ZOMHz  scope  (model  4224)  in  combination  with  a  Dell  Vostro  Laptop. 
Data  is  acquired  at  10  MSamples/s  with  a  precision  of  12  bits.  The  PicoScope  and  laptop  create  a 
raw  data  file  which  is  less  than  SOOkB  in  size,  by  detecting  and  storing  only  periods  in  which 
there  are  pulses  (ie  energy)  in  the  data.  The  system  and  file  format  are  documented  in  [5].  From 
2009  to  2010  the  testing  was  done  using  a  National  Instruments  PXI-5122  digitizer  and  PXI  1022 
platform,  sampling  at  2  MSamples/s  with  14  bit  precision.  Because  the  entire  data  stream  was 
recorded,  this  produced  data  files  in  excess  of  27MB  making  memory  management  a  critical 
consideration. 


1 .4  Characteristics  of  the  T aser  X26  Waveform 


This  section  provides  details  on  the  waveform,  definitions  and  specifications  of  the  parameters 
of  interest  for  a  TASER  X26.  The  X26  pulse  consists  of  an  "Arc  phase"  and  "Main  phase"  as 
shown  in  fig.  2.  The  pulses  are  delivered  in  a  cycle  consisting  of  approximately  95  pulses  over  5 
seconds,  at  the  rate  of  approximatelyl9  pulses  per  second.  In  a  X26  pulse  waveform,  the 
information  which  electrically  characterises  the  weapon  is  primarily  derived  from  the  Main 
phase,  where  most  of  the  pulse  energy  resides.  The  Main  phase  delivers  approximately  100  pC 
of  charge,  whereas  the  arc  phase  has  only  about  10  pC. 


Voltage 

or 

Current 


Figure  2  Taser  X26  pulse  waveform  showing:  (i)  Peak  voltage  (or  peak  current)  computed  using 
the  waveform  Main  phase,  (ii)  Net  charge  calculated  as  the  absolute  value  of  area  under  the 
current  curve  in  the  Main  phase  (y),  (Hi)  Full  pulse  duration  starting  at  an  initial  crossing  of  ±50V 
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and  finishing  when  pulse  decreases  to  ±50V,  (iv)  Monophasic  Charge  calculated  as  the  maximum 
of  absolute  values  of  y  and  d,  (v)  Total  charge  computed  as  the  sum  of  the  absolute  values  of  y 
and  &. 

Table  1  shows  the  parameters  of  individual  pulses  which  were  identified  in  [3].  The  list  of  pulse 
parameters  was  partially  based  on  the  definitions  by  Taser  International  [6] 


Table  1  Definition  of  Parameters 


Parameter 

Definition 

Peak  voltage/current  using  the  Main 
phase: 

The  peak  voltage  or  peak  current  of  an  X26  pulse  is 
computed  using  the  Main  phase  of  a  voltage  or  current 

curve. 

Net  charge  using  the  Main  phase: 

The  net  charge  is  computed  using  the  Main  phase  of 
an  X26  pulse.  The  Net  charge  is  calculated  as  the  area 
under  the  Main  phase  section  starting  at  the  first 
crossing  of  the  Main  phase  above  OV^  (OviA)  and 
finishing  when  pulse  decreases  to 

83mA  =  50F/600a 

Pulse  duration  using  full  pulse: 

Pulse  duration  of  an  X26  pulse  starting  at  an  initial 
crossing  of  ±50K  and  finishing  when  pulse  decreases 
to  +50F. 

Pulse  repetition  rate: 

The  pulse  repetition  rate  is  computed  by  1)  counting 
the  number  of  pulses  in  a  full  cycle,  2)  measuring  the 
cycle  time  from  the  start  of  the  first  pulse  to  the  start 
of  the  last  pulse  and  3)  the  pulse  repetition  rate  is  one 
less  than  the  total  number  of  pulses,  divided  by  the 
cycle  time. 

Monophasic  Charge: 

The  Monophasic  Charge  is  calculated  as  the  maximum 
of  the  absolute  values  of  y  and  0,  where  y  is  the 
integral  of  the  positive  current  in  the  X26  pulse 
waveform  and  0  is  the  integral  of  the  negative  current 
in  the  X26  waveform,  respectively 

Total  Charge: 

The  total  charge  is  calculated  as  the  sum  of  the 
absolute  values  of  y  and  0. 
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2  Carleton  System 


2.1  Test  Environment  at  Carleton  University 

The  Carleton  University  CEW  test  system  was  originally  designed  in  2009  and  has  seen 
incremental  improvements  over  the  past  two  years.  Its  design  was  based  on  the  following 
guidelines: 

1.  Testing  equipment  is  to  be  as  mobile  and  portable  as  possible  to  allow  testing  to  take  place 
in  front  of  the  client  and  to  minimize  out-of-service  time  for  the  CEW. 

2.  At  the  same  time,  the  test  environment  should  preserve  the  integrity  of  measurement  and 
analysis  by  adhering  to  sound  engineering  practices.  This  means  that  the  components  of  the 
test  system  must  tolerate  movement  in  secure  Pelican-type  cases.  The  equipment  and 
system  as  a  whole  must  be  configured  and  calibrated  in  the  field  as  easily  as  it  would  in  a 
static  laboratory  environment.  Other  variables  such  as  temperature  and  humidity  were 
measured  and  noted  in  each  separate  environment. 

The  test  environment  for  in-situ  testing  is  always  indoors  in  a  heated  and/or  air-conditioned 
secure  space.  Temperatures  are  generally  20C  ±  2C,  and  humidity  is  generally  between  20%  and 
40%.  All  test  spaces  are  controlled  areas  where  there  is  no  public  access  and  "technical 
tourism"  is  held  to  a  minimum,  and  thus  physical  distractions  are  minimized. 

Weapons  are  controlled  by  police  officers  and  handed  to  the  Testing  Engineer  by  the  officer  for 
insertion  in  the  mount.  The  CEW's  are  returned  to  their  cases  or  storage  mounts  when  not 
being  physically  examined  or  electrically  tested.  There  are  no  special  requirements  for 
ventilation  or  air  handling  in  the  test  environment  because  there  are  no  harmful  odours  or 
particulate  emissions  for  the  firing  of  the  weapon.  There  is  no  projectile  from  firing  and  there  is 
no  noise  from  firing  into  a  calibrated  load  so  there  are  no  special  requirements  with  respect  to 
position  of  firing  officer  or  safety  equipment  such  as  goggles  or  ear  protection.  The  software  and 
hardware  used  in  the  testing  environment  are  controlled  by  the  test  engineer.  There  is  no 
unauthorized  use  or  modification  of  either  hardware  or  software  outside  the  purpose  of  testing 
CEW's  for  which  it  was  intended. 

2.2  Carleton  University  CEW  analysis  software 

CEW  testing  at  Carleton  University  is  performed  in  two  logically  separate  steps,  as  outlined  in  fig 

3.  First  data  are  acquired  using  an  acquisition  hardware  system.  Data  are  digitized  and  stored  to 
a  file  (the  file  format  is  documented  in  the  next  section,  2.3).  Next,  the  data  are  analysed  using 
the  analysis  software  documented  here.  CEW  pulse  data  from  the  raw  data  file  are  loaded,  and 
pulse  parameters  calculated. 
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Figure  3  Schematic  diagram  of  workflow  ofCEW  testing  at  Carleton.  Step  1:  the  CEW  is 
connected  to  the  electrical  test  system  and  data  from  a  weapon  firing  Is  captured  and  stored  to  a 
raw  datafile.  Step  2:  the  raw  datafile  Is  loaded  by  the  CEW_analyse  software  and  parameters 
are  calculated. 

2.3  Carleton  University  CEW  output  file  format 

Data  acquired  by  the  CEW  data  capture  hardware  are  written  to  a  file  with  a  format 
documented  in  [2]  and  in  Table  2.  This  file  format  is  designed  to  save  storage  space  by  storing 
only  pulse  data  during  the  time  that  the  weapon  is  active.  Between  pulses,  when  the  output 
voltage  is  zero,  the  signal  value  is  not  stored. 

Table  2  Raw  Data  file  format 


Data  Contained 

Data  Type 

Valid  Data 

Header  data 

File  Format 

String 

CU-CEWdata  VER  1.00 

Date  and  time  of 

Data  collection 

String 

2010/10/26:14h34 

Serial  Number 

String 

X00181701 

Description 

String 

A  string  of  length  < 
65535  bytes,  (length, 
string) 

Owner's  name 

String 

A  string  of  length  < 
65535  bytes,  (length, 
string) 

Comments 

String 

A  string  of  length  < 
65535  bytes,  (length, 
string) 
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Gain 

Float  (4  bytes) 

Positive  real  values 

Offset  (V) 

Float  (4  bytes) 

Positive  real  values 

Resistance  (Q) 

Float  (4  bytes) 

Positive  real  values 

Sampling  Rate  (FIz) 

Float  (4  bytes) 

Positive  real  values 

Number  of  Pulses 

Integer  32-bit 

0  to  2^^  -  1 

Pulse  length 

Integer  32-bit 

0  to  2^^  -  1 

Interpulse  Time 
(From  the  beginning 
of  a  pulse  to  the 
beginning  of  the  next 
pulse) 

ti(s) 

Integer  16-bit 

Oto  2^®  -  1 

t2(s) 

Integer  16-bit 

0  to  2^^  -  1 

Integer  16-bit 

Oto  2^®  -  1 

tn-i(s) 

Integer  16-bit 

0  to  2^^  -  1 

Interpulse  Time 
(From  the  end  of  a 
pulse  to  the  end  of 
the  next  pulse) 

ti(s) 

Integer  16-bit 

0  to  2^^  -  1 

t2(s) 

Integer  16-bit 

0  to  2^^  -  1 

Integer  16-bit 

0  to  2^^  -  1 

tn-l(s) 

Integer  16-bit 

0  to  2^^  -  1 

Sampling  Data 

Si,i(V) 

Integer  16-bit 

-2^^  to  2^^  -  1 

Integer  16-bit 

-2^^  to  2^^  -  1 

S  Pulse  length,  1  (V) 

Integer  16-bit 

-2^^  to  2^^  -  1 

Sl,2(V) 

Integer  16-bit 

-2^^  to  2^^  -  1 

Integer  16-bit 

-2^®  to  2^^  -  1 

S  Pulse  length,  2  (V) 

Integer  16-bit 

-2^^  to  2^^  -  1 

Sl.n(V) 

Integer  16-bit 

-2^^  to  2^^  -  1 

Integer  16-bit 

-2^®  to  2^^  -  1 

S  Pulse  length,  n  (V) 

Integer  16-bit 

-2^^  to  2^^  -  1 

2.4  Carleton  University  CEW  anaiysis  software 

The  software  to  perform  the  analysis  of  CEW  waveforms  is  called  CEW_analyse.m.  It  is  designed 
to  be  run  under  Matlab  or  GNU  Octave  interpreters,  and  is  used  to  analyse  signals  from  CEW's 
acquired  on  the  PicoScope.  It  comprises  382  lines  and  13  functions.  Table  3  sets  out  the  names 
of  the  functions  and  what  each  function  does,  as  well  as  some  ancillary  information  such  as  line 
numbers,  variables  passed  in,  other  functions  called  by  the  function  being  described.  The 
MatLab  code  can  be  seen  in  annex  A. 
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Table  3  CEW_analyse.m  function  descriptions 


Name  of  function 

Lines 

What  it  does 

Calls  other  function: 

1 

CEW_analyse 

12-65 

For  each  data  file  entered  the  data  is 
read  and  placed  in  to  data  and  ss,  finds 
the  params  for  each  data  point  plots 
the  voltage  at  each  point  if  requested 
and  then  call  print_outputs 

print_outputs 

read_data 

calibrate_data 

retrigger_data 

find_pts 

find_patarms 

print_output 

2 

print_outputs 

70-173 

Outputs  the  first  two  tables  on  page  13 
of  the  test  procedure,  for  each  shot 
calculates  interpulse  length,  total  firing 
length,  interpulse  distance,  median 

PRF,  PRF  from  last  second,  set  flag  if 

PRF  <16.5  prints  a  string  header  of 
table.  Prints  the  third  table  on  page  13 
of  the  test  procedure. 

max_m  i  n_a  vg_a  vg8 

outlines 

prt_log 

3 

outlines 

177-193 

Define  column  widths  and  string 
headers  for  tabular  output  in  FITML 
(a.k.a  page  13  output) 

4 

find_pts 

197-220 

Defines  3  points  on  the  waveform: 
start,  zcross,  endpoint 

my_medfilt 

5 

find_patarms 

225-242 

Defines  limits  of  arc  phase  and  main 
phase  and  produces  dT.  Uses  Vpk  to 
calculate  Ipk,  PD,  Qnet,  Qmono,  Qtotal 

6 

system_settings 

244-262 

System  settings  for  Picoscope 

7 

getstr 

265-267 

Reads  name  from  the  files  and  gets  a 
string  of  that  length 

8 

read_data 

271-305 

Reads  text  data  from  data  file.  Reads  9 
system  settings  variables  from  data  file. 

If  no  pulses  detected  print  warning.  If 
max_data  exceeds  clip  limits,  print 
warning.  Print  string  identifying  which 
data  is  being  accessed 

get_str 

prt_log 

9 

calibrate_data 

311-321 

Calibrates  the  data  to  the  gains 
associated  with  the  equipment  being 
used  and  calculates  the  voltage  from 
the  data.  Determines  the  time  at  each 
sampling  point. 

10 

retrigger_data 

326-344 

Reduces  the  length  of  data  from  250  ps 
to  200  ps 

11 

prtjog 

347-354 

Creates  an  output  file  and  prints  to  it 

12 

my_medfilt 

358-365 

Applies  a  median  filter,  each  location 
the  data  is  replaced  with  the  median 
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value  of  the  data  points  next  to  it. 

13 

max-min-avg-avgS 

370-382 

Finds  absolute  maxima,  absolute 
minima,  mean  and  average  of  last  8 
pulses.  Sets  flag  if  last  8  pulses  exceed 
max  value  or  are  less  than  min  value. 
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2.5  Physical  schema  of  the  software  data  elements 


Table  4  describes  the  variables  used  by  each  function  in  CEW_analyse.m.  The  description  of 
each  variable  is  given  in  Annex  B. 

Table  4  Physical  schema  overview 


Function 

Variables  passed  in 

Variables  at  Output 

Local  Variables 

CEW_analyse 

f_name, options 

ss,  data,  V,  t,  pts, 
params 

print_outputs 

params,ss 

interpulse,  PRF, 
mmaa,  flag 

outlines 

0 

find_pts 

Vi,ss 

pts 

start,  zcross,  endpt 

find_patarms 

Vi,  t,  pts,  ss 

PP 

arc_phase, 
main_phase,  dT 

system_settings 

ss 

getstr 

fid 

str 

len 

read_data 

f_name  ,ss 

data,  ss 

fid 

calibrate_data 

datajnput,  ss 

data,  V,  t 

retrigger_data 

V,  t,  ss 

Vo,  t 

thresh,  Vi, 

prtjog 

str,  varargin 

txt,  resultsfile 

my_medfilt 

s,  n 

s 

Is,  idx 

max-min-avg-avgS 

data,  minval,  maxval 

m_m_a_a8,flag 

2.6  Logical  Data  Model 

This  section  describes  the  logical  data  model  of  the  software  and  its  data  elements  of  the 
CEW_analysis.m  software.  Figure  4  gives  the  underlying  logical  data  model  with  which  we  were 
working  for  the  design  of  the  software. 
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Figure  4  Schematic  representation  of  the  Logical  Data  Model. 

2.7  Business  rules  incorporated  into  the  software 

2.7.1  Business  Rules 

The  business  rules  used  by  Carleton  in  the  analysis  of  a  CEW  are  as  follows: 
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♦  Equipment  is  setup  and  tested  at  the  facility  in  which  testing  is  to  be  conducted. 

♦  The  weapon  owners  are  briefed  about  the  test  procedure  to  be  followed. 

♦  The  Carleton  Pico-scope  acquisition  software  is  executed  on  the  PC. 

♦  The  Description  and  Owner  of  the  CEWs  being  tested  is  entered 

♦  The  individual  testing  of  each  individual  CEW  begins  and  is  repeated  for  each 
weapon 

■  The  serial  number  and  comments  for  the  CEW  being  tested  is  entered 

■  Comments  include  any  visual  observations  or  special  equipment  such  as 
handgrip  camera 

■  The  program  waits  for  the  firing  to  begin 

■  The  trigger  is  pulled 

■  The  recording  begins  when  the  voltage  crosses  lOOmV;  this  voltage  is  after  a 
1000:1  voltage  reducing  probe 

■  System  records  each  pulse  from  the  weapon,  starting  at  50ps  before  this 
voltage  crossing  the  trigger  (using  the  pre-trigger  feature  of  the  hardware)and 
200ps  after  the  voltage  trigger. 

■  A  trigger  level  of  crossing  (either  +50V  or  -50V)  is  used  for  all  pulses 

■  Then  the  time  between  each  pulse  is  measured 

■  Recording  of  pulses  stops  if  either  of  the  following  conditions  are  met:  a)  120 
pulses  are  recorded,  or  b)  the  testing  user  hits  a  key.  In  almost  all  cases, 
recording  is  stopped  by  the  user. 

■  The  program  then  saves  data  to  a  raw  data  file  in  the  format  outlined  in 
section  2.3 

♦  The  operator  ensures  that  the  raw  data  file  and  analysis  code  are  placed  in  the 
same  folder. 

♦  The  analysis  code  is  then  called  in  Matlab  using  the  following  call; 

CEW_analyse({file  name  one;  file  name  two;  file  name  three},  options).  File  names 
are  with  extensions  and  in  single  quotations,  e.g.  'binlSout- 

(PoliceDepartment,X00-123456,97,21)-01276110295cvt.dat' 

♦  For  each  pulse  in  the  data  stream,  the  parameters  are  calculated,  as  defined  by. 
Data  points  refers  to  individual  sampled  voltage  values  (after  correcting  for  the 
gain  and  offset): 

■  Peak  voltage:  determined  by  the  maximum  data  point  on  the  data  points  in  the 
Main  phase  (note:  this  is  not  the  maximum  absolute  value). 

■  Peak  current:  determined  by  the  peak  voltage  divided  by  the  load  resistance. 

■  Pulse  duration  is  the  time  from  the  pulse  start  point  to  the  pulse  end  point.  It  is 
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calculated  by  subtracting  (pulse  end  point)  -  (pulse  start  point),  in  which  the 
points  are  calculated  as: 

♦  Pulse  start  point:  the  time  at  which  the  signal  crosses  from  V>-  50  V 
to  V<-50  V. 

♦  Pulse  end  point:  the  last  time  at  which  the  filtered  signal  crosses 
from  V>50V  to  V<50V.  Filtering  is  a  median  filter  in  which  each  point  is 
the  median  of  the  samples  within  a  width  of  2  us  in  each  direction. 
The  last  crossing  is  identified  by,  first,  identifying  all  crossings  of  the 
50V  threshold  by  the  filtered  signal  in  each  pulse,  and,  then,  taking  the 
last  one. 

■  Main  Phase:  The  main  phase  is  the  time  from  the  main  phase  start  to  the  pulse 
end  point  (defined  above). 

♦  The  main  phase  start  is  the  first  time  the  signal  crosses  from  V<0V  to 
V>0V  after  the  pulse  start  point 

■  Net  Charge  is  calculated  by  summing  the  data  points  in  the  main  phase 
multiplied  by  the  sampling  period  and  divided  by  the  resistance  (note:  this  is 
the  sum,  not  the  sum  of  absolute  values) 

■  Monophasic  charge  is  calculated  by  the  maximum  of  A  and  B,  where:  A  is  the 
sum  of  the  absolute  value  of  all  data  points  with  V>0  during  the  pulse  duration 
multiplied  by  the  sampling  period,  and  divided  by  the  resistance;  and  B  is  the 
sum  of  the  absolute  value  of  all  data  points  with  V<0  during  the  pulse  duration 
multiplied  by  the  sampling  period,  and  divided  by  the  resistance 

■  Total  charge  is  calculated  by  summing  the  absolute  value  of  data  points  during 
the  pulse  duration  multiplied  by  the  sampling  period,  and  divided  by  the 
resistance. 

■  Interpulse  time  for  pulse  N  is  calculated  as  the  time  from  pulse  N  start  point  to 
pulse  N+1  start  point. 

♦  Based  on  the  parameters  calculated  for  each  pulse,  summary  parameters  are 

calculated  for  the  entire  pulse  train  for  the  weapon  firing. 

■  For  parameters  Peak  voltage.  Peak  current,  Net  Charge,  and  Pulse  Duration, 
the  following  summary  parameters  are  calculated 

♦  Maximum,  the  maximum  of  all  pulse  values 

♦  Minimum,  the  minimum  of  all  pulse  values 

♦  Average,,  the  average  (mean)  of  all  pulse  values 

♦  Average-TI,  the  average  (mean)  of  the  last  8  pulse  values 

■  For  parameters  Total  Charge,  and  Monophasic  Charge,  the  following  summary 
parameters  are  calculated 

♦  Maximum,  the  maximum  of  all  pulse  values 
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♦  Minimum,  the  minimum  of  all  pulse  values 

♦  Average,  the  average  (mean)  of  all  pulse  values 


■  Based  on  the  parameter  interpulse  time,  the  a  pulse  rate  is  calculated  for  each 
pulse  as  l/(interpulse  time).  Using  this  value,  the  following  summary 
parameters  are  calculated 

♦  Maximum,  the  maximum  of  all  pulse  values 

♦  Minimum,  the  minimum  of  all  pulse  values 

♦  Average,  the  average  (mean)  of  all  pulse  values 

♦  Burst  Length,  the  sum  of  all  interpulse  time  values 

♦  Average-TI,  the  (Number  of  pulses  -  1)  /  (burst /engfh) 

♦  Summary  parameters  are  compared  to  the  specified  values  in  the  appendix  A 
(TASER  M26)  and  appendix  B  (TASER  X26)  of  [1].  Values  outside  of  the  specified 
range  are  identified  (with  the  string  '****') 

♦  Individual  pulse  and  summary  parameters  are  saved  to  a  file 

♦  Summary  parameters  are  displayed  to  the  user  and  saved  to  a  file 


2.8  Additional  data  elements  and  business  rules  that  will  be 
required  to  fully  implement  the  data  requirements  of  CEW 
Test  procedure 

The  current  version  of  CEW_analyse.m  and  test  apparatus  used  by  Carleton  in  the 
testing  of  CEW  requires  some  additional  data  elements  to  conform  the  Carleton  test 
procedure,  this  can  be  seen  in  Table  5. 

5.  Requirements  to  conform  to  the  Carleton  test  Procedure 


Data  Element 

Comments 

Temperature 

It  is  currently  being  recorded  by  hand  on  to  paper.  It  should  be  in  the 
raw  data  file  with  the  other  system  settings 

Humidity 

It  is  currently  being  recorded  by  hand  on  to  paper.  It  should  be  in  the 
raw  data  file  with  the  other  system  settings 

Remaining  battery 
capacity 

It  is  not  currently  being  recorded;  it  has  been  seen  to  have  no  effect  on 
how  the  weapon  performs.  It  should  be  recorded  and  placed  into  the 
raw  data  file  with  the  other  system  settings  so  as  to  conform  to  the 

DRDCCSSCR  2011-18 


23 


Carleton  test  procedure. 

Software  Version 

It  is  not  currently  being  recorded.  It  should  be  recorded  to  conform  to 
the  Carleton  test  procedure.  In  practice,  most  police  departments  are 
buying  the  latest  DPM  which  contains  the  most  recent  software  and 
cycling  this  DPM  through  the  inventory  of  CEW's. 

Owner  of  the  CEW 

The  owner  of  the  CEW  is  currently  in  the  raw  data  file  but  is  not  being 
consistently  recorded  while  testing. 
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Summary  and  Conclusions 


This  report  describes  the  test  analysis  software  developed  at  Carleton  University  (CEW)  that  was 
written  to  support  the  work  on  the  Conducted  Energy  Weapons  (CEW)  Test  Procedure,  [3].  The 
objectives  of  the  report  is  to  document  the  software,  its  structure,  business  rules,  and  to  identify 
any  data  elements  required  to  fully  support  the  published  test  procedure.  Within  the  report,  each 
of  these  objectives  is  systematically  described. 

This  report  is  designed  to  support  the  Conducted  Energy  Weapons  Strategic  Initiative  (CEWSI), 
which  plans  to  test  a  number  of  CEWs  and  do  analysis  on  legacy  test  data  with  the  objective  of 
better  understanding  and  interpreting  the  performance  of  the  devices,  with  a  goal  of  developing 
an  updated  CEW  Test  Procedure  and  providing  advice  to  Canadian  Eaw  Enforcement  on  these 
devices. 
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Annex  A  MATLAB  code  for  CEW analyse.m 


1  %%This  code  is  designed  to  analyse  CEW  X26  weapons  according  to  the 
test 

2  %%procedure  found  online  at : 

3  %%http : / / curve . car let on . ca/papers/ 2  010/ CEW-Test-Procedure-2 010- 
verl . 1 . pdf 

4  %%(c)  A.  ADLER  and  0.  MARSH,  2009-2011 

5  %%This  code  may  be  copied  under  the  GPL  licence  version  2  or  3  found 
at : 

6  %%http : / / WWW . gnu . org/li censes/ gpl. html 

7 

8  %%Main  function  of  CEW_analyse  reads  data  from  the  file  in  f_name 
and  places 

9  %%it  in  data  and  ss,  checks  for  the  presence  of  energy  in  the  file, 

10  %%finds  V  at  all  sample  points,  plots  all  voltage  values  and  repeats 
for  the 

11  %%files  entered  in  to  the  function  then  calls  print_outputs 

12  function  CEW_analyse { f_name,  options); 


13 


14  ss  =  system_settings ; 

15 

16  if  nargin>l ; 

17  ss. plotting  =  options . plotting; 

18  ss . debug_plot  =  options . debug_plot ; 

19  ss . use_f liter  =  options . use_f liter ; 

20  ss.open  =  options . open; 

21  else 

22  ss. plotting  =  1; 

23  ss . debug_plot  =  0; 

24  ss . use_f liter  =  0; 

25  ss.open  =  1; 

2  6  end 

27 

28  if  ss. plotting;  clf;  end 

29 

30  n_shot  =  size { f_name) ; 

31  n_shot  =  n_shot{l,l); 

32  if  n_shot>3; 

33  fprintf{'Too  many  files  entered  only  first  three  will  be  usedXn'); 

34  n_shot  =  3 

35 

36  elseif  n_shot  ==  2; 

37  ss=  [ss,  ss]; 

38  elseif  n_shot  ==  3; 

39  ss=  [ss,  ss,  ss]; 

4  0  end 

41 

42  for  j  =  l:n_shot; 

43 

44  if  ss { j ). plotting;  figure (j);  end 

45 

46  [data,  ss{j)]  =  read_data {  char { f_name { j ) ) ,  ss{j)); 


47 


48 


[data,V,t]  =  calibrate_data {  data,  ss{j)); 
[V,t]  =  retrigger_data (V, t  ,ss{j)); 


49 

50 

51 


for  i=l:ss{j) .pulses 
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52 

Vi  =  V( : , 1) ; 

53 

if  all(Vi==0);  break;  end 

54 

pts  =  find_pts (  Vi,  ss(j)  ); 

55 

params (j,  i)  =  find_patarms ( 

Vi,  t. 

pts,  ss  (  j)  ) ; 

56 

if  ss(j) .plotting 

57 

11=  [pts . start , pts . zcross , 

pts . endpt ] ; 

58 

plot (t , Vi , t ( 11 ) ,  Vi(ii),'*' 

) ; 

59 

if  i==l;  hold  on; end 

60 

end 

61 

end 

62 

end 

63 

print_outputs (  params, ss) 

64 

fprintf ( ' Calculations  completeXn ' ) ; 

65 

if  ss(l) .open;  open ( ' summary_results . html 

'  ) ; end 

66 

67 

68 

%%print_outputs  prints  the  results 

of  the 

calculations  in  html 

format 

69  %%which  follows  the  format  of  page  13  in  CEW  test  procedure 

70  function  print_outputs (  params,  ss) 

71  tablel2  =  strcat { ' <table  width  =  "800"  border="2"  rules  =  "all" 
bordercolor  =" 

black"><tr> ' ,  .  .  . 

72  ' <td  width  =  "500"><b>  Conductive  Energy  Weapon  Test  Report  <b/>  <td 
width  =  "300"> ' ,  .  .  . 

73  'Date:  %s  <tr>  <td  >  Weapon:  <td  >  Serial  Number:  %s<tr>  <td  > 

Police 

Service :  %s  '  ,  .  .  . 

74  '  <td  >  Police  Officer:  <tr><table/><br> ' ,  .  .  . 

75  ' <table  width  =  "800"  border="2"  rules  =  "all"  bordercolor  ="black" 
><tr> ' , . . . 

76  '<td  width  =  "200"  >Visual  Inspection<td  width  =  "600">  Case  Battery 
Electrodes ' , . . . 

77  '<tr><td>  Data  Download  Performed  <td><tr><td>  Comments 
<td>%s<tr> '  ,  .  .  . 

78  ' <td>  Software  Version  <td><tr><td>  Battery  Charge  <td><tr>',  ... 

79  ' <td>  Battery  Model  and  Serial<td><tr><td>  Temperature  <td><tr><td> 
Humidity<td> ' , . . . 

80  '<tr><td>  Atmospheric  Pressure  <td><table/><br> ' ) ; 

81 

82  prt_log {tablel2 ,  ss{l) .meastime, ss (1) . serial , ss { 1 ) .owner, 
ss { 1 ) . comment ) ; 

83  clear  tablel2 

84 

85 

86  n_shot  =  size  (params) ; 

87  n_shot  =  n_shot{l); 

88  interpulse  =  zeros (size  (params)); 

89  for  i=l:n_shot 

90  interpulse ( i , 1 : length ( ss ( 1 ). interpse_e ( 2 : end-1 )) )  = 
ss(i) . interpse_e ( 2 : 

end-1 ) ; 

91  f iringlength_2=  sum ( interpulse ( 1 ,  :  )  )  ; 

92  firinglength  =  f iringlength_2 *( ss ( 1 ). pulses )/( ss ( 1 ). pulses-2 )  + 

mean 

( [params (1,  : )  .pulse_duration] )  ; 

93 

94  medPRF=  1 /median ( interpulse ( 1 ))  ; 

95  PRF  =  (ss(i) .pulses-1)  /  firinglength; 

96 

97  if  PRF<16.5;  flag  =  '****'; 
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98  else  flag  =  '  '  ;  end 

99 

100  prt_log { ' Shot  %d:  %d  pulses  in  %4.2f  s.  Pulse  Rep  Freq  =  %4.2f  p/s 

"O 

s<br><br> ' ,  ... 

101  i,  ss  { 1 ). pulses ,  f iringlength,  PRF,  flag); 

102 

103  end 

104 

105  oo  =  outlines; 

106 

107  %creating  Third  output  table 

108  tables  =  strcat { ' <table  width  =  "800"  border="2"  rules  =  "all" 
bordercolor  =" 

black"><tr> ' , . . . 

109  '<td  width  =  "150"><th  width  =  "150"  colspan=" %d" >  Max  <th  width  = 
"150"  colspan="%d"> '  ,  . . . 

110  'Min<th  width  =  "150"  colspan=" %d" >  Avg<th  width  =  "150"  colspan="% 
d">  Avg  -  TI '  ,  .  .  . 

111  '<td  width  50>  SPEC  <tr><td  >  Firing  No'); 

112 

113  tables  =  sprintf (tableS,  n_shot,  n_shot,  n_shot,  n_shot) ; 

114 

115  for  1=1:4 

116  for  j  =  l:n_shot 

117  tables  =  sprintf (strcat (tableS,  '<td  align  =  "center">%d' ) , j) ; 

118  end 

119  end 

120 

121  tables  =  strcat (tableS ,  '<td>'); 

122 

123  prt_log (tables ) ; 

124 

125  flag.pv  =  flag. pc  =  flag.nc  =  flag.pd  =  flag.pr  =''; 

126  flag.tc  =  flag. me  = 

127  for  i  =  l:n_shot 

128  [mmaa(i, :) .pv,tflag]  =  max_min_avg_avg8 ( 

[params (1,  : )  . peak_voltage ]  , 

1400,2520  ) ; 

129  flag.pv  =  strcat ( flag . pv,  tflag) ; 

130  [mmaa(i,  :)  .pc,  tflag]  =  max_min_avg_avg8 ( 

[params ( 1 ,:). peak_current ]  ,  2.3, 

4.2  )  ; 

131  flag. pc  =  strcat ( flag . pc,  tflag); 

132  [mmaa(i,  :)  .nc,  tflag]  =  max_min_avg_avg8 ( 

[params (1,  : )  .net_charge] *le6,  80, 

125  )  ; 

133  flag.nc  =  strcat ( flag . nc,  tflag); 

134  [mmaa(i, :) .pd, tflag]  =  max_min_avg_avg8 ( 

[params (1, :) .pulse_duration] 

*le6, 105,  155  ) ; 

135  flag.pd  =  strcat ( flag . pd,  tflag); 

136  [mmaa(i, :) .pr  , tflag]  =  max_min_avg_avg8 ( 

[ 1 . /nonzeros (interpulse  ( 1 ,  :  ) )  ] , 

0,  inf  ) ; 

137  flag.pr  =  strcat ( flag . pr,  tflag); 

138  [mmaa(i,  :)  .tc,  tflag]  =  max_min_avg_avg8 ( 

[params (1, :) .total_charge] *le6, 

0,  inf  ) ; 

139  flag.tc  =  strcat ( flag . tc,  tflag); 
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140  [mmaa{i,  :)  .mc,tflag]  =  max_min_avg_avg8 { 

[params{i, :) . monophasic_charge ] 

*le6,  0,  180  )  ; 

141  flag. me  =  streat { f lag . me,  tflag) ; 

142  end 

143 

144  prt_log (oo .pv)  ; 

145  prt_log (oo . table,  mmaa.pv); 

146  prt_log (oo . f lag,  flag.pv); 

147 

148  prt_log (oo .pe)  ; 

149  prt_log (oo . table,  mmaa.pe); 

150  prt_log (oo . f lag,  flag.pe); 

151 

152  prt_log (oo . ne)  ; 

153  prt_log (oo . table,  mmaa.ne); 

154  prt_log (oo . f lag,  flag.ne); 

155 

156  prt_log (oo .pd)  ; 

157  prt_log (oo . table,  mmaa.pd); 

158  prt_log (oo . f lag,  flag.pd); 

159 

160  prt_log (oo .pr)  ; 

161  prt_log (oo . table,  mmaa.pr); 

162  prt_log (oo . f lag,  flag.pr); 

163 

164  prt_log (oo . te)  ; 

165  prt_log (oo . table,  mmaa.te); 

166  prt_log (oo . f lag,  flag.te); 

167 

168  prt_log (oo .me)  ; 

169  prt_log (oo . table,  mmaa.me); 

170  prt_log (oo . f lag,  flag. me); 

171 

172  prt_log ( ' </table> ' )  ; 

173  return 

174 

175 

176  %%outlines  defines  the  formatting  of  the  tabular  HTML  output 

177  funetion  o  =  outlines; 

178  o. table  =  ' <td  align  =  "eenter">  %.2f'; 

179  o.flag  =  '<td  align  =  "eenter">  %s ' ; 

180  %  peak  voltage  (main  phase) 

181  o.pv='<tr><td  >  Peak  Voltage  (V) '; 

182  %  peak  eurrent  (main  phase) 

183  o.pe=  '<tr><td  >Peak  Current  (A)  '; 

184  %  net  eharge  (main  phase) 

185  o.ne=  '<tr><td  >Net  Charge  (&#181;C) '; 

186  %  pulse  duration  (full  pulse), 

187  o.pd=  '<tr><td  >Pulse  Duration  (&#181;s) '; 

188  %  pulse  repitition  rate 

189  o.pr=  '<tr><td  >Pulse  Rep  Rate  (Hz) '; 

190  %  Monophasie  Charge 

191  o.me=  '<tr><td  >Mono .  Charge  (&#181;C) '; 

192  %  Total  Charge 

193  o.tc=  '<tr><td  >Total  Charge  (&#181;C)  '; 

194 

195  %%find_pts  finds  the  location  of  three  points  on  the  waveform,  the 
start 

196  %%point,  first  zero  crossing,  and  the  end  point 

197  function  pts  =  find_pts (  Vi,  ss  ); 
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198  Vi  =  my_medfilt {  Vi,  0); 

199  start  =  find {Vi<=-50) ; 

200  start  =  start (1); 

201 
202 

203  zcross=  f ind {Vi>=0 ) ; 

204  zcross { zcross<start ) =  [  ] ; 

205  zcross  =  zcross  (1); 

206 

207  Vi  =  my_medfilt {  Vi,  20); 

208  endpt  =  find{Vi  <  50); 

209  delay  =  30e-6*ss . sample_f ;  %  30  us 

210  endpt (endpt <  (zcross  +  delay))  =  []; 

211  if  isempty (endpt) 

212  endpt  =  190e-6*ss . sample_f-l 

213  else 

214  endpt=  endpt (1); 

215  end 

216 

217 

218  pts. start  =  start-1; 

219  pts.zcross=  zcross; 

220  pts. endpt  =  endpt; 

221 

222  %%f ind_patarms  uses  the  data  for  each  pulse  to  calculate  peak 
voltage, 

223  %%peak  current,  net  charge,  pulse  duration,  monophasic  charge  and 
total 

224  %%charge 

225  function  pp  =  find_patarms (  Vi,  t,  pts,  ss  ); 

226  arc_phase  =  Vi (pts . start :  (pts. zcross  -1)); 

227  main_phase  =  Vi (pts . zcross : pts . endpt ) ; 

228  dT  =  mean (dif f (t) ) ; 

229 

230  %  peak  voltage  (main  phase) 

231  pp . peak_voltage  =  max (main_phase) ; 

232  %  peak  current  (main  phase) 

233  pp . peak_current  =  pp . peak_voltage  /  ss  .  resistance; 

234  %  net  charge  (main  phase) 

235  pp . net_charge  =  sum(  main_phase  )  *  dT  /  ss . resistance; 

236  %  pulse  duration  (full  pulse), 

237  pp . pulse_duration  =  t (pts. endpt)  -  t (pts . start ) ; 

238  %  Monophasic  Charge 

239  pp . monophasic_charge  =  sum(  Vi  .*  (Vi>0)  ) *dT  /  ss . resistance; 

240  %  Total  Charge 

241  pp . total_charge  =  sum(  abs (  Vi  )  ) *dT  /  ss . resistance; 

242 

243  %%system_settings  defines  the  default  settings  of  the  testing 
apparatus 

244  function  ss=  system_settings 

245  ss . sample_f=10e6; 

246  ss  .  resistance  =  595; 

247  ss.gain  =  1/1000; 

248  ss. offset  =  0;  %actually  -6mV 

249  ss .probegain=1000; 

250  ss . input_gain=l  ; 

251  ss . clip_limit  =  2^(12-1)  -  1; 

252  ss . f ilef ormat  = 

253 

254  ss.meastime  = 

255  ss . serial  =  ' ' ; 
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I  }  . 


256  ss . description= 

257  ss. owner  =  ' ' ; 

258  ss. comment  =  ' ' 

259  ss. pulses  =  0; 

260  ss.pulse_len  =  0; 

261  ss . interpse_s  =  0; 

262  ss . interpse_e  =  0; 

263 

264  %%getstr  reads  a  length  of  a  string  for  a  file  and  then  the  string 

265  function  [str]  =  getstr(fid) 

266  len=  f read { f id, 1 , ' int 1 6 ' ) ; 

267  str=  char (f read (fid, [l,len], ' char*l ' ) ) ; 

268 

269  %%read_data  reads  the  data  and  system  settings  from  the  file  and 
checks 

270  %%if  there  is  data  in  the  file  and  if  signal  is  vertically  clipped 

271  function  [data,ss]  =  read_data {  f_name,  ss) ; 

272  fid  =  fopen { f_name, ' rb ' ) ; 

273  ss . file format  =  char (f read (fid, [1,20], ' char* 1 ' ) ) ; 

274  if  ss  .  fileformat  (end)  ~=0;  error  {' format '),’ end 

275  ss.meastime  =  getstr { f id) ; 

276  ss. serial  =  getstr { fid)  ; 

277  ss . description=  getstr { fid) ; 

278  ss. owner  =  getstr { fid)  ; 

279  ss. comment  =  getstr { fid) ; 

280  ss.gain  =  f read { fid, 1 ,' float ') ; 

281  ss. offset  =  f read ( fid, 1 ,' float ')  ; 

282  ss . resistance  =  f read { fid, 1 ,' float ') ; 

283  ss.sample_f  =  f read { fid, 1 ,' float ') ; 

284  ss. pulses  =  f read { fid, 1 ,' int32 ') ; 

285  ss.pulse_len  =  f read { fid,  1 ,' int32 ')  ; 

286  ss . interpse_s  =  f read ( fid, ss . pulses ,' float ') ; 

287  ss . interpse_e  =  f read { fid, ss . pulses ,' float ') ; 

288  data  =  f read { fid,  [ ss . pulse_len, ss . pulses ],' int 1 6 ') ; 

289 

290  f close (fid) ; 

291  if  ss.pulse_len  ==500 

292  ss . clip_limit  =  32764; 

293  end 

294  % - Analysing  the  stripped  data - 1- 

295  max_data  =  max (data); 

296  min_data  =  min (data); 

297  if  max (abs { [max_data, min_data] )) <1000 

298  prt_log {' WARNING :  No  pulses  detected  in  data  file<br>'); 

2  99  data=  [  ] ; 

300  return 

301  elseif  {max_data>=ss . clip_limit )  | {min_data<-ss . clip_limit ) 

302  prt_log { ' <br>WARNING :  Signal  is  vertically  clipped<br> ' ) ; 

303  end 

304 

305  prt_log { ' <br>Data  loaded:  f ile=%s<br><br> ' , f_name) ; 

306 

307 

308  %%calibrate_data  calibrates  the  data  to  the  gains  associated  with 
the 

309  %%testing  apparatus  and  calculates  the  voltage  and  time  at  each 

310  %%sampling  point 

311  function  [data,V,t]  =  calibrate_data {  data_input,  ss) ; 

312  if  ss . use_f liter 

313  [data,  ip]  =  func_despike_phasespace3d {  double {data_input) ,  1,2 

)  ; 
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314  else 

315  data  =  double {data_input) ; 

316  end 

317 

318 

319  data=  { ss . input_gain) *data; 

320  V  =  {data*ss . gain+ss . of fset) *ss .probegain; 

321  t  =  { 0 : ss . pulse_len-l ) /ss . sample_f ; 

322 

323  %%retrigger_data  takes  the  250us's  of  sample  points  and  reduces  to 
200us 

324  %%of  data  sample  points.  It  takes  lOus  before  the  pulse  begins  and 
1 90us 

325  %%after  it  begins 

326  function  [Vo,t]  =  retrigger_data (V,  t,  ss); 

327  thresh  =  100; 

328  lowlim  =  intl6 {10e-6*ss . sample_f ) ; 

329  VN  =  size (V, 2) ; 

330  idx  =  -lowlim: 190e-6*ss . sample_f-l; %spanning  200us 

331  Vo  =  zeros { length ( idx) , VN) ; 

332 

333  for  i=l : VN 

334  Vi  =  V{ : , i) ; 

335  ff  =  find{  abs{Vi)  >  thresh); 

336  if  length { ff) >=1 

337  ff=  ff  (1) ; 

338  Vo{:,i)  =  Vi {  ff+idx  ); 

339  lim  =  intl6 {lle-6*ss . sample_f ) ; 

340  Vo { lowlim : lim,  i )  =  linspace{Vo{lowlim,i),Vo{lim,i),lim- 
(lowlim-1) ) ; 

341  end 

342  end 

343 

344  t=  t  { 1 : length { idx) )  ; 

345 

346  %%prt_log  outputs  the  string  called  in  to  the  function  to  output 
file 

347  function  prt_log{str,  varargin  ); 

348  txt=  sprintf (str, varargin) :}) ; 

349 

350 

351  resultsfile  =  f open { ' summary_results . html '  ,  ' at ' ) ; 

352  fprintf (resultsfile,  txt); 

353  fclose (resultsfile)  ; 

354 

355  %%my_medfilt  is  a  median  filter  to  array  s  entered  into  the  function 

356  %%each  data  point  is  replaced  with  the  median  data  point  of  the  2n+l 

357  %%data  points  to  either  side  of  it 

358  function  s  =  my_medfilt (  s,  n) ; 

359  Is  =  length(s); 

360  idx=  ones (2*n+l,l)*(l:ls)  +  (-n:n)' *ones (1,1s) ; 

361  idx(idx<l)  =  1; 

362  idx(idx>ls)  =  Is; 

363  s  =  reshape ( s ( idx) ,  [],  Is) ; 

364  s  =  median (s,l); 

365 

366  %%max_min_avg_avg_avg8  takes  data  and  determines  the  absolute 
minimum, 

367  %%absolute  maximum  and  mean  of  all  pulses.  Then  determines  mean  of 
the 
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368 

the 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 
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%%last  eight  pulses.  It  will  throw  a  flag  up  if  it  is  greater  than 
%%maximum  or  minimum  values . 

function  [m_m_a_a8 , f lag]  =  max_min_avg_avg8 {  data,  minval,  maxval  ) 
m_m_a_a8{l)  =  max (data); 
m_m_a_a8{2)  =  min (data); 
m_m_a_a8{3)  =  mean (data); 
m_m_a_a8{4)  =  mean (data (end- { 0 : 7 ) ) ) ; 

flag  =  ' ' ; 

i f  m_m_a_a8 { 4 ) >maxval ; 

f lag=  ' *  '  ; 
end 

i f  m_m_a_a8 { 4 ) <minval ; 

f lag=  ' # ' ; 
end 
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Annex  B  Physical  schema  of  the  software  data 
elements 


Data  Elements  used  by  CEW_analyse: 


Data 

Element 

What  it  is 

Type 

Example 

f_name 

The  File  name  of  the 

file  to  be  analysed 

A  array  of  strings 

CEW-X00526878- 

(Thu_Dec_16_15.27.47_2010).dat 

options 

An  Object  contain  the 
options  for  it  plotting, 
debup_plot,  use  filter, 
and  open  results  file 

An  object  with  4 

Boolean  values 

contained  in  it. 

options. plotting  =  0 
options.debug_plot  =  0 
options.use_filter  =  0 
options.open  =  0 

ss 

The  meta  data  for  the 

testing 

An  object  with  17 
fields  with  variety 

of  values  stored 

ss.meastime  = 

'2010/10/26:14h34'; 
ss.serial  =  'x00526878'; 
ss.description=  'ops'; 

data 

Raw  data  from  the  file 

An  array  of  int-16 

values 

assorted  values  ranging  from 
-2^^  to  2^^  -  1 

V 

Data  calibrated  to  the 

associated  gains 

An  array  of  int-16 

values 

assorted  values  ranging  from 
-2^^  to  2^^  -  1 

t 

Times  at  each  point  in 

the  wave  form 

An  array  of  floating 
point  numbers 

values  ranging  from  Os  to  200us 
spaced  evenly 

pts 

The  points  on  the 

waveform  at  with  the 

wave  starts,  crosses 

the  zero  point  for  the 

first  time  and  the  end 

point 

An  object 

containing  3  integer 

values. 

Values  range  from  0  to  2000 

params 

An  array  containing 

the  characteristics  for 

each  waveform  in  the 

pulse  train 

An  array  of  objects 
with  each  object 
containing  6  float 

fields 

params(l).peak_voltage  =  1850 
params(l).peak_current  =  3.1 
params(l).net_charge  =  105.3e-6 
params(l).pulse_duration  = 

125.2e-6 

params(l).monophasic_charge  = 

110.4e-6 
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params(l).total_charge  =  IZO.le- 

6 

Data  elements  used  by  print_outputs:  params  and  ss  are  equivalent  to  the  fields  in 
CEW_analyse. 


Data  Element 

What  it  is 

Type 

Example 

interpulse 

Values  of  the  time 

between  pulses 

An  array  of  float 

values 

values  ranging  from  0.04s 

to  0.2s 

PRF 

the  Pulse  repetition 
frequency  of  the  firing 

An  float  value 

accepted  values  range  from 
16.5-20pps 

mmaa 

values  of  the  mean, 
median,  average  and 
average  of  the  last  8  for 
each  parameter 

An  array  of  an  object 
containing  7  fields 
each  being  an  array 

of  floats 

mmaa.pv  =  2000.3  1600.4 

1850.5  1840.7 

Flag 

flag  for  if  any  parameter 

is  out  of  tolerance 

An  object  with  7 

fields  with  are 

strings 

flag.pv  =  '****' 
flag.nc  =  '####' 

Data  elements  used  by  outlines: 


Data  Element 

What  it  is 

Type 

Example 

0 

The  definition  for  the 

html  output  for  each 
parameter 

An  object  containing 

9  fields  each  being  a 
string. 

o.table  =  '<td  align  = 

"center">  %.2f' 

o.pv  =  '<tr><td>Peak 

Voltage  (V)' 

Data  elements  used  by  find_pts:  ss  and  pts  are  equivalent  to  the  fields  in  CEW_analyse. 


Data  Element 

What  it  is 

Type 

Example 

Vi 

The  voltage  data  for 
the  pulse 

An  array  of 
integers 

Assorted  values  ranging 
from  -2^^to2^^-l 

start 

The  index  in  the  array  at 
which  the  pulse  starts 

An  integer 

Normal  values  range  from 

0-1000 

zcross 

The  index  in  the  array  at 
which  the  pulse  crosses 

zero  for  the  first  time  in 

the  waveform 

An  integer 

Normal  values  range  from 

0-1000 
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endpt 

The  index  in  the  array  at 

An  integer 

Normal  values  range  from 

with  the  pulse  ends 

1000-2000 

Data  elements  used  by  find_patarms:  ss,  pts  and  t  are  equivalent  to  the  fields  in  CEW_analyse. 
Vi  is  equivalent  to  the  field  in  find_pts. 


Data  Element 

What  it  is 

Type 

Example 

PP 

The  parameters  (pv, 
pc,  nc,  pd,  me  and  tc) 
of  the  waveform 

An  object  with  6 
float  fields. 

pp.peak_voltage  =  1850 
pp.peak_current  =  3.1 
pp.net_charge  =  105.3e-6 
pp.pulse_duration  =  125. 2e- 
6 

pp.monophasic_charge  = 
110.4e-6 

pp.total_charge  =  120. le- 
6 

arc_phase 

the  voltage  values  for 
the  data  points  in  the 
arc  phase  of  the  pulse 

An  array  of  int-16 

Assorted  values  ranging 
from  -2^^  to  2^®  -  1 

main_phase 

the  voltage  values  for 
the  data  points  the  main 
phase  of  the  pulse 

An  array  of  integers 

Assorted  values  ranging 
from  -2^®  to  2^®  -  1 

dT 

The  time  increments  for 

between  the  data  points 

A  float  value 

O.lus 

Data  elements  used  by  system_settings:  ss  is  equivalent  to  the  field  in  CEW_analyse 


Data  elements  used  by  getstr: 


Data  Element 

What  it  is 

Type 

Example 

fid 

File  id  of  the  file  being 
read  from 

File  id 

N/a 

str 

The  string  taken  from 

the  file 

A  string 

'X00526878' 

len 

Length  of  the  string 
being  read  from  the  file 

A  int-16 

9 
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Data  elements  used  by  read_data:  f_name,  ssand  data  are  equivalent  the  field  in  CEW_analyse. 
fid  is  equivalent  the  field  in  getstr. 


Data  elements  used  by  calibrate_data:  ss,  data,  V  and  t  are  equivalent  the  field  in  CEW_analyse. 


Data 

Element 

What  it  is 

Type 

Example 

datajnput 

The  raw  data  which  needs 

to  be  calibrated 

An  array  of  int- 
16  values 

assorted  values  ranging 
from  -2^^to2^^-l 

Data  elements  used  by  retrigger_data:  V,  t  and  ss  equivalent  the  field  in  CEW_analyse.  Vi 
equivalent  to  the  field  in  find_pts. 


Data  Element 

What  it  is 

Type 

Example 

Vo 

The  data  points  after  the 

data  has  been 

retriggered 

An  array  of  int-16 

Assorted  values  ranging 
from  -2^^  to  2^®-l 

thresh 

The  threshold  voltage  to 

determine  when  the 

pulse  begins 

An  int-16  value 

100 

Data  elements  used  by  prtjog: 


Data  Element 

What  it  is 

Type 

Example 

str 

The  string  to  be 
outputted 

A  string 

'<td  align  =  "center">%d' 

varargin 

Any  values  required  to 
be  outputted 

Any  number  of 
values  of  any  type 

1 

txt 

the  output  string  with 

the  values  added 

A  string 

'<td  align  =  "center">T 

resultsfile 

the  file  id  of  the  output 

file 

file  id 

N/A 

Data  elements  used  by  my_medfilt: 


Data  Element 

What  it  is 

Type 

Example 

s 

The  data  points  which 

An  array  of  int-16 

Assorted  values  ranging 

38 


DRDCCSSCR  2011-18 


need  to  be  filter 

from  -2^^  to  2^^-l 

n 

Number  of  points  on 

each  side  to  take  the 

median  from 

An  int-16 

The  normal  value  is  40 

Is 

Number  of  data  points  in 

s 

An  int-16 

Normally  2000 

Data  elements  used  by  max-min-avg-avgS: 


Data  Element 

What  it  is 

Type 

Example 

data 

The  parameter  values  for 
each  pulse 

An  array  of  floats 

Assorted  values 

minval 

The  lower  spec  limit 

A  float 

For  Peak  Voltage  it  is  1400 

maxval 

The  upper  spec  limit 

A  float 

For  Peak  Voltage  it  is  2520 

m_m_a_a8 

The  mean,  median, 
average  and  average-8 

values 

An  array  of  floats 

Assorted  values 

flag 

Flag  if  any  parameters 
are  out  of  spec 

A  string 

>(c  ^ 
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Annex  C  Algorithms  for  Each  Function 
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5.  find_patarms(Vi,t,pts,ss)  lines  225-242 


define  arcphase 
and  mainphase 


6.  system_settings  lines  244-262 


7.  getstr  lines  265-267 


Hard  eode  magie  numbers  for  Pieoseope 
and  new  load 
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12.  prt_log(str,varargin)  lines  347-354 


Creates  an  output  file  and  prints  to  it 
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14.  max_min_avg_avg8(data,  minval,  maxval)  lines  370-382 


Finds  absolute  minima,  maxima, 
mean  and  average  of  last  8 
pulses 


Sets  flag  if  last  8  pulses  exeeed 
max  value  or  are  less  than  min 
value. 
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List  of  symbols/abbreviations/acronyms/initialisms 


CEW 

Conducted  Energy  Weapon 

CEWSI 

Conducted  Energy  Weapon  Strategic  Initiative 

CSS 

Center  for  Security  Sciences 

DRDC 

Defence  Research  &  Development  Canada 

TI 

TASER  International 
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